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Abstract 

Heat shock response is a stress response to temperature changes and a consecutive increase 
in amounts of unfolded proteins. To restore homeostasis, cells upregulate chaperones facilitating 
protein folding by means of transcription factors (TF). We here investigate two heat shock systems: 
one characteristic to gram negative bacteria, mediated by transcriptional activator cr 32 in E. coli, 
and another characteristic to gram positive bacteria, mediated by transcriptional repressor HrcA 
in L. lactis. We construct simple mathematical model of the two systems focusing on the negative 
feedbacks, where free chaperons suppress cr 32 activation in the former, while they activate HrcA 
repression in the latter. We demonstrate that both systems, in spite of the difference at the TF 
regulation level, are capable of showing very similar heat shock dynamics. We find that differences 
in regulation impose distinct constrains on chaperone-TF binding affinities: the binding constant 
of free cr 32 to chaperon DnaK, known to be in 100 nM range, set the lower limit of amount of free 
chaperon that the system can sense the change at the heat shock, while the binding affinity of HrcA 
to chaperon GroE set the upper limit and have to be rather large extending into the micromolar 
range. 
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Introduction 



Cellular homeostasis is essential for proper protein folding and function. The perturba- 
tions to homeostasis, e.g. due to change in temperature or osmotic pressure, result in protein 
unfolding or/and misfolding. Heat shock, i.e., sudden increase of temperature, causes such 
protein unfolding and misfolding and can result in cell death. To counteract the heat shock, 
cells upregulate production of chaperons and proteases - enzymes that help folding the un- 
folded proteins and degrade terminally misfolded proteins, respectively. The heat shock 
response is one of the stress responses characteristic for nearly all living organisms. Interest- 
ingly, the protein sequence of most chaperones and proteases is well conserved from bacteria 
to humans It is, however, unclear if the features of heat shock response are also preserved 
at the level of the architecture of regulatory circuits governing heat shock response. 

In this article we attempt to answer this question and derive useful insights by comparing 
the heat shock in E. coli and L. lactis These organisms utilize two different mechanisms; a 
system with a 32 and DnaK in E. coli and a system with HrcA and GroE in L. lactis. Both 
mechanisms are widely observed in microorganisms. A transcriptional activator RpoH, a 32 
homolog, is found in the alpha-, beta-, and gamma-proteobacteria, while a transcriptional 
inhibitor, HrcA, is widely distributed in eubacteria but not in the gamma-proteobacteria. In- 
terestingly, there also exist bacteria which have both systems and, furthermore, a regulatory 
loop between a 32 and HrcA is predicted in some beta-proteobacteria [2J. 

Heat shock responses have been extensively studied both experimentally and theoretically. 
In experiments, protein sequences and their regulatory mechanisms are revealed in the both 
systems (figure [T]) [3|. While a 32 system is modeled to quite large extent 0, to our 
knowledge there is no modeling work on HrcA system. Our aim in this study is to construct a 
simple model based on known experimental data for each system and theoretically investigate 
similarities and differences in the regulatory features and the dynamical responses mediated 
by a 32 and HrcA. 

One of the striking similarities is emerging at the level of the dynamics of the transcription 
regulators: a 32 and HrcA. Both systems respond with a sharp peak in the rate of production 
of new chaperones: upon a temperature shift, a fast increase up to 4-5 fold within 5-10 
minutes (corresponding to about 0.1 generation time [5]) is followed by a rapid decline to 
a new steady state that is about 1.5 fold of the one at the starting temperature in both 
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FIG. 1: Comparison of the reaction mechanisms between (A,C) cr 32 -DnaK system and (B,D) HrcA- 
GroE system. Views for gene regulations and protein reactions (A and B) and outline illustrations 
for distinct reactions (C and D). 

a 32 and HrcA system 5|, |6|- This similarity is particularly interesting as the mechanisms of 
transcriptional regulation are very different: while a 32 is an activator, HrcA is a repressor 
(see figure [I]) . 



Model 



In this section, we explain how we construct our models based on existing experimental 
observations. 



a -DnaK system 

Our model with a 32 and DnaK is in large similar to that outlined in [4] with exception 
of the differences outlined below. The main players in the model are: a 32 (er; transcription 
factor), DnaK (D; chaperon), and unfolded proteins (U). (In the following we will denote 
total concentration of protein X as [X f ] and the concentration of free protein as [X *].) 

a 32 is unstable and is present only in small amounts ([a 1 ] ~ 200 proteins per cell). 
Under steady state conditions it is sequestered by chaperones, such that the most a 32 exists 
as a complex ([a*] ~ [o"-D])- For simplicity, we assume it is produced at a constant rate 
independent of temperature. In addition, for a simpler comparison with HrcA system, we 
do not include the stabilization of a 32 (half-life changes from 1 to 8 min) during the heat 
shock. These two are the main differences from the model in . 

When not bound to DnaK, a 32 forms a complex with RNA polymerase (RNAp) and 
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targets RNAp for transcription of heat shock proteins, including DnaK. For simplicity we 
will refer to this [RNAp : a 32 } complex as "free" <r 32 , [a*\ (i.e. not bound to DnaK). 

Being a chaperone, DnaK facilitates proteins folding and thus forms transient complexes 
with unfolded proteins. A temperature shift destabilizes the folding of existing folded pro- 
teins and also hinders folding of de novo synthesized proteins, thus sequestering all the 
existing chaperones and creating the demand for additional chaperones. The demand for 
additional chaperones is sensed and mediated by cr 32 : as long as there are enough unfolded 
proteins to keep chaperones sequestered away from a 32 , it will facilitate transcription of heat 
shock proteins. Namely, the regulatory network has a negative feedback loop (see figure [TA 
and C), i.e., <r 32 activates DnaK by transcriptional activation ("slow" reaction) while DnaK 
inhibits a 32 by complex formation ("fast" reaction). 

We assume the reactions a' + D* ^ oD and U* + D* ^ UD are in equilibrium as 
the kinetics of complex formation and dissociation between cr 32 and DnaK and between 
unfolded proteins and DnaK are much faster than transcription and translation. Based on 
these observations and assumptions, we describe the system's dynamics with the following 
equations; 

[^]= V - ls [a t ]- lc [aD], (1) 
[Ut] = F{T)- lus [UD], (3) 

m = vm, ( 4) 

with the conditions of mass conservation 

[a 1 ] = [a f ] + \oD] ~ 200 nM, (6) 
[£>*] = [D f ] + [uD] + [UD] - 20, 000 nM, (7) 
[U t ] = [Uf] + [UD]. (8) 

The rate of change in a 32 (eq. (QQ)) is given by constitutive transcription rate, rj, dilu- 
tion due to cell division with the cell doubling, l/j s , and degradation with fast rate, j c . 
Since a 32 is degraded mainly through chaperone-dependent degradation by FisH [7j, the 



fast degradation term depends on the complex [<jD]. The production of chaperon DnaK 
(first term in eq. fl2])) transcriptionally activated by free a 32 , and parameterized by 
the maximal production rate a^- For [a?] dependence, we adopt the Michaelis-Menten form 
with dissociation constant K a . 

Chaperones are typically stable proteins with the half life comparable to cell doubling 
time 1/js, which is included in the second term of eq. (E}. The first term in eq. ([3]), F(T), 
represents that the production of unfolded protein is increasing function of the temperature, 
T. Unfolded proteins can re-fold correctly with help of chaperons and are thus removed 
at a rate proportional to [UD] with the rate 7„ s (the second term in eq. (J3J)). Eqs. (jl]) 
and (E]) represents complex formations in equilibrium, with Kj and Ky. are the dissociation 
constants between free a 32 and free DnaK and between free unfolded protein and free DnaK, 
respectively [16 j. 



We fix some of the parameter values according to the experimental observations as follows. 
We set 7 C to a unit time (= 1), which shows fast degradation of a complex aD as most a 32 
exist as the complex in the steady state. The time scale of the fast degradation, l/j c is 
assumed to be around 1 min js] [l7|. l/j s is a time scale for slow degradation of a 1 and D* 
and it is set to the inverse number of cell division time (~ 30). We estimate rj ~ 200 nM/min 
from an observation that ler'l ~ 200 nM. The dissociation constant between chaperones and 
a 32 range between 5 /iM and 19 nM between correspondingly DnaK/DnaJ and a 32 |7|. We 
have set the constant to be Kj = 100 nM, following the choice in (4]. 

The rest of parameters, aca, K a , Kk,j U s, and F(T), have been chosen so that the model 
reproduces experimental observations, i.e. (i) [D*] ~ 20,000 nM [s,[lo|]. (ii) DnaK produc- 
tion rate changes 4 ~ 6 times for a peak and its new steady state after heat shock becomes 

fi n 

1.5 times of the before heat shock [5[. (iii) The peak time is less than 5 time units [18( and 
the peak shape of the DnaK production rate is symmetric (iv) The steady state 

amount of free unfolded proteins should be kept small both before and after heat shock. This 
is affected by Kh,j us , and F(T), the properties related to unfolded proteins. In this paper 
we fix Kk = 1 nM so that Uf is in nanomolar range in steady state, and fit the other two 
parameters. We tested higher values of (up to = 1,000 nM, which would bring Uf 
to micromolar range) and they all give a proper response as long as F(T) and 7 US adjusted 
properly to account for the timing of the peak. 
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HrcA-GroE system 



Next, we construct a model for HrcA and GroE system, where GroE (G) is a chaperon 
and HrcA (H) transcriptionally represses GroE. We adopt the reaction mechanism shown 
in figure 2 in llj ; HrcA repressor is released from the ribosomes as an inactive protein (H l ), 
which can not bind to its operator, and it has to interact with the GroE chaperonin system 
to become active (H a ). The inactive HrcA (H l ) interacts with chaperon GroE {HG) and 
the active HrcA is released. The active HrcA (H a ) is able to bind to its operators and 
transcriptionally inhibits the production of GroE, while the active HrcA becomes inactive 
again at a constant rate, i.e., upon dissociation from its binding site, HrcA is in its inactive 
form again [12|. In this model, we assume that the total amount of HrcA ([H l ] = [H a ] + 
[H % ] + [HG]) is a constant for simplicity. As in case with DnaK, GroE chaperons makes a 
complex with an unfolded protein and helps it to refold correctly. Similar to er 32 system, we 
assume that the two reactions (H % + G* ^ HG and U* + G* ^ UG) are fast compared 
with other reactions and always in the equilibrium states. 

This system also has a negative feedback loop between transcription factor HrcA and 
chaperon GroE. However, the regulation is opposite; the active HrcA inhibits GroE with 
a slow reaction (transcriptional inhibition) and GroE activates the inactive HrcA with a 
fast reaction (enzymatic modification). From the points described, we obtain the following 



reaction equations; 



[H']= /y/fG]- 7c [/P], (9) 

^=^TTwm h - lAG% (10) 

[U']=F{T)- luh [UG], (11) 

[HG] = M, (12, 



K, 



K, 

with the conditions of mass conservation 



= [H a ] + [H l ] + [HG], (14) 
[G f ] = [G f ] + [HG] + [UG], (15) 
[U'] = [Uf] + [UG]. (16) 
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Eq. (Q, represents the time evolution of [H a ] with constant total amount of HrcA, where 
the production rate of [H a ] is given by activation from H l through forming a complex 
[HG] with the rate @h, and inactivation happens with a fast rate j c . The time evolutions of 
chaperon [G*] (eq. (I10p) and of unfolded protein [U*] (eq. (jlip ) are similar to eqs. (|2J) and (J3J), 
except that [G 1 ] is transcriptionally inhibited by [H a ] in eq. (I10B . Here, g is the maximum 
production rate of GroE, Kh is the dissociation constant of active HrcA to GroE promoter, 
and '-fuh is the rate of removal of unfolded proteins by GroE. Eqs. (fT2]) and ([TBI represents 
complex formations in equilibrium, with K\ and K m are the dissociation constants between 
inactive HrcA and free GroE and between free unfolded protein and free GroE, respectively 

Most of the parameters were not experimentally measured for the HrcA-GroE system. 
However, for a fair comparison between the two systems, we, whenever possible, use same 
parameter values as in a 32 -DnaK system (see Tabled for correspondence relation). Thus we 
assume similar concentrations for transcription factor [H f ] = 200 nM (~ [a*]) and chaperones 
[G 1 ] ~ 20,000 nM (~ [D 1 ]). Although these numbers have not been validated experimen- 
tally, they match the typical concentrations of transcription factors (of order 100 nM) and 
chaperones (of order 10 /xM). 



The fastest doubling time in L. lactis is also about 30 min [13], thus we use the same 
doubling time as in E. coli with a corresponding doubling time 7 S = 0.03. Unlike a 32 , HrcA 
is a stable protein (in B. subtilis), with the half-life more than 60 minutes [l4]. However, it 
has been suggested that HrcA is present in two conformations, one is active and another is 
inactive, and the equilibrium between these two states is modulated by GroEL, which shift 
the equilibrium towards the active state 12] . Thus 7 C , which we set to be 1, is representing 
the rate of conversion from active to inactive state, rather than protein half-life as is the 
case in cr 32 . 

The time of the peak, fold induction at the peak and fold change of the new steady state 
are overall similar in the activity of CIRCA operon [6J (corresponding to GroEL production 
rate) and er 32 governed chaperone (DnaK) production. This allows us to use the same 
criteria (i) - (iv) to set K m , (3 g , Kh,j u h, Ph, F(T), and Ki. We study the response of this 
model assuming the same conditions as for a 32 system (figure [2] C,D). (In HrcA model, we 
assumed [H l ] is constant for simplicity. However, the model also works even if we discard 
this postulation and include that HrcA inhibits its own transcription 11], i.e., the time 
evolution of [H l ] will be described as [#*] = /i/(l + \H a ]/K t ) - %[&].) 
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cr 32 -DnaK system 


HrcA-GroE system 


V 


200 


[H*] 


200 


7c 


1 




1 


7s 


0.03 




0.03 


K k * 


1 




1 


K i 

J 


100 


Ki * 


100,000 


ad * 


2,000 


Pg * 


2,000 


K s * 


10 


K h * 


10 


7ns * 


0.5 


luh * 


0.5 






Ph * 


3 


F(T) * 


3000 -> 9000 


* 


3000 -»• 9000 



TABLE I: Parameter values used in the model. * marks parameters chosen such as to reproduce 
rapid transient response measured in j], Qj.The rest of the parameters are based on experimental 
data as described in the text. 



Results 



Model predicts much weaker binding affinity between HrcA and chaperones. 



In figure [2] we show that both <r 32 -DnaK and HrcA-GroE systems are able to reproduce 
experimental observations. Production of chaperones shows characteristic sharp peak with 
a fast increase up to 4 — 6 fold within 5 minutes and a following decline to a new steady 
state that is about 1.5 fold of the pre-stimulus one. When choosing unknown parameters, 
our initial strategy was to use the same values for corresponding parameters in each of the 
systems (see Tabled]). Remarkably, this was possible for all but one parameter: the binding 
affinity of TF to chaperons, Kj and K\. It appears that while a 32 binds tightly to the 
chaperons (Kj = 100 nM), it is essential that HrcA is bound only weakly with a micromolar 
binding constant (K\ = 100 fiM). 

We next study a response upon an inverse heat shock, i.e., a response when temperature 
is suddenly decreased. Inverse heat shock response has been studied experimentally in E. 



coli 



15, 



15] (to our knowledge no data exist for L. lactis) and is characterized by a rapid 



decrease in chaperone production with a consequent slow increase to a new steady state 
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FIG. 2: Heat shock response in cr 32 -DnaK system (A,B) and HrcA-GroE system (C,D). Heat 
shock is induced at t = 0. (A,C) shows the time evolution of the chaperone production rate 
normalized by the pre-stimulus level. (B,D) shows the time evolution of the density of total 
chaperones(dotted lines), total unfolded-protein(bold lines), and complex formed by chaperon and 
unfolded proteins (broken lines). Parameters used in simulations are shown in Tabled 



that is lower than before temperature decrease. 

In figure [21 we show the chaperon production rate upon an inverse heat shock. The 
response is simulated using the parameter values in Table U except for F(T) which is reversed 
(suddenly decreased at t — 0.). Both models for cx 32 -DnaK and HrcA-GroE systems showed 
very similar responses, which fit well with experimental results; the chaperon production 
rate shows a rapid transient decrease and recovers slower compared with a direct heat shock 
response. The fact that the model works without specific tuning of parameters to the inverse 
heat shock supports that our simple models hit the essential points of the actual reaction 
mechanisms. 

Why the recovery to the new steady state is slower in the inverse heat shock in both 
systems? The explanation naturally emerges from our model: as chaperons are stable pro- 
teins, the only way to recover to a new steady state upon decrease in unfolded proteins 
is by dilution due to cell division. Thus this slow time scale for the recovery is given by 
l/j a or the time scale of cell division in both systems. In case of direct heat shock the 
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time to recover to new steady state (right after the peak, once there are enough chaperons 
produced) is governed by the turnover time for active TFs, 7 C , which is much faster than 
the rate governed by cell doubling time j s . 
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FIG. 3: Inverse heat shock response in <r 32 -DnaK system (A,B) and HrcA-GroE system (C,D). 
(A,C) shows the time evolution of chaperone production rate normalized by the pre-stimulus level. 
(B,D) shows the time evolution of the density of total chaperones (dotted lines), total unfolded- 
protein(bold lines), and complex formed by chaperon and unfolded proteins (broken lines). Param- 
eters used in simulations are shown in Table U except for production rate of unfolded proteins, 
F(T), which decreases at t = 0, F(T) = 9, 000 -)■ 3, 000. 



The difference in negative feedback architectures requires different constraints on 
TF-chaperon binding affinities 

The reaction mechanisms of the two systems resemble each other in that there exists a 
negative feedback loop between a transcription factor and a chaperon. However, loops are 
organized such that TF is an activator in one and is an inhibitor in another. In the following 
we will demonstrate how this difference leads to the distinctly different binding affinities of 
TF to chaperons. 

The constrains on binding affinities can be understood when we look at how TF and 
chaperons are related in steady state. From eqs. ([T]), (jlj) and (jSJ) we obtain the expression 
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for free a to be 



a 



V K j _ v/ls 



lsKj + ( lc + ls )[Df] i+ m 

7s-Kj/(7c+7«) 



(17) 



This is a decreasing function of the free chaperon [D*], and [er*] approaches a constant value 
77/7^ when [D^] < j s Kj/(^ c + 7 S ). Here, j s Kj/(^ c + 7 S ) ~ 3 as we fix to 7 S = 0.03 and 
Kj = 100 based on experimental observations. 

Similarly for the HrcA, from eq. fl9]), f fT2|) and f fT4|) . we find 

[ J 7c^ + (7c + 4)[G'/] 1+ -|L • 1 > 

This is an increasing function of the free chaperon [G*], and [H a ] approaches a constant 
value (3h[H t ]/(j c + (3h) when [Gf] ^> r y c Ki/{^ c + (3h). Note that both K\ and are unknown 
parameters for the HrcA system. 

These expressions (with parameters from Table [T|) are plotted in figure HJA. as functions 
of free chaperons. As chaperons activate HrcA and inhibit a 32 , HrcA is increasing and a 32 
is decreasing with increasing amounts of free chaperons. Each of the curves has two char- 
acteristic regimes: a) Insensitive regime where TF are insensitive to changes in chaperone 
concentration. This corresponds to a nearly flat region in the plot, where concentration of 
TF does not depend or depends very weakly on the chaperone concentration (chaperones< 3 
for a 32 , and chaperones> 10 5 for HrcA) and b) Sensitive regime where change in chaper- 
one concentration results in a change in TF concentration (chaperones> 3 for a 32 , and 
chaperones< 10 5 for HrcA). 

For the system to be responsive and adjust the production rate of chaperons (controlled by 
amounts of active TFs) in response to changes in unfolded proteins (reflected by the amounts 
of free chaperons), it is essential for the system to function within sensitive regimes. The 
peculiar feature of the two systems is that the insensitive regimes lie in the opposite ends of 
chaperone concentrations; a 32 system is sensitive as long as the free chaperon \D'\ is above 
the threshold concentration, while HrcA system can work as long as the free chaperon [G*] 
is below the threshold concentration. At the same time, the maximal concentration of free 
chaperon is limited by the amount of total proteins (~ 20, 000), and the minimum represents 
the case when all chaperones are bound to unfolded proteins. As in steady state the amounts 
of free chaperones vary between 5, 000-10, 000, it is crucial that the sensitive regime spans 
this range. As insensitive regime for a 32 lies in the range of small concentrations (threshold 
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jsKj/^jc + 7 S ) ~ 3 with experimentally evaluated parameters), it will always be in the 
sensitive range. On the contrary, in HrcA system, the threshold must be larger than the 
typical steady state concentrations of free chaperons, i.e. r y c Ki/( , y c + 0h) > 10,000. This 
impose the condition on binding affinity to be large enough, Ki > 10, 000, because the 
threshold value mainly depends on K\ and tuning other parameters such as 0h or [H 1 ] does 
not affect threshold values much (figure HJ3 and C). 

It now becomes clear why we could not obtain the responses in HrcA with K\ ~ Kj ~ 100, 
as this tight binding in HrcA system would decrease the sensitivity regime to be below 100 
nM or below typical steady state concentrations of free chaperons. On another hand, for a 
similar reason, we can not use weak binding affinity for a 32 , i.e. Kj = 100, 000, as it will 
shrink the sensitivity region to be above 3, 000 nM. 

Alternatively, if in reality K\ ~ 100 nM this would imply that the sensitivity regime is 
very narrow, which means that the steady states of free chaperones have to vary between 1 
and 100 nM. In principle this could be the case, however, this would imply that the system 
is not robust to sudden increases in unfolded proteins. 

This difference in the chaperone-TF dissociation constant between the two systems is 
very critical. While both systems have negative feedback as a core regulatory mechanism, 
our model predicts that the details of how each of the feedbacks is realized result in very 
different dissociation constants. One can test this prediction experimentally by varying bind- 
ing affinity in HrcA system or sensitivity domains and investigate how this affects response 
dynamics. 

Parameter Robustness 

The results presented so far were based on a single set of parameters, chosen to reproduce 
experimental data. Some of them are fixed to a known experimental values as already 
mentioned in the model section, while the rest of the parameters are fitted to reproduce 
heat shock dynamics; there are 5 fitting parameters for cr 32 -DnaK system and 7 for HrcA- 
GroE (see Table [T|). To understand how constrained is our parameter choice we studied 
robustness of parameters. Figure [5] shows how much a given parameter can be changed (the 
maximal fold of change from the values shown in Table H]) with preserving proper heat shock 
response. 
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FIG. 4: Relation between free chaperon and transcription factor in a steady state given by eq. (|17p 
and eq. (|18p . [af] is with the parameters in Table [J (bold broken line), while for [H a ] some lines 
changing a parameter [Ki in A, [H t ] in B, and j3h in C) is shown and fih = 3, [H t ] = 200, and 
K\ = 100, 000 are used if not otherwise specified 



We choose a proper response to the one characterized by (A) a peak with more than 2 and 
less than 10 fold change (normalized to the pre-stimulus level) (B) occurring within 10 unit 
times and (C) recovering to a new steady state that is less than 2 fold of the pre-stimulus 
one in chaperon production rate. Interestingly, most parameters can only be changed at 
most a few fold for both systems. 

One of the main reasons why system is rather sensitive to parameter choice is due to 
the importance of the stoichiometry between chaperons and unfolded proteins: if either 
chaperons or unfolded proteins are in excess there will be no peak. Excess of chaperons (e.g. 
due to high chaperone production rate g ) will absorb a sudden increase in unfolded 
proteins, so that there will be no increase in chaperone production and thus no peak. On 
the other hand excess of unfolded proteins (e.g. low ad, (3 g ) will result in a state where 
chaperone production is maximally activated already before the shock. Thus a further 
increase in unfolded proteins will not lead to the increase in production rate of chaperons. 
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FIG. 5: Parameter robustness measured as fold change over the reference parameter set (shown in 
Table|I]) that is able to reproduce a reasonable heat shock response (defined by conditions (i)-(iv)). 
E.g. the value of 1 on y axis correspond to the values shown in Table HI Red marks results for 
cr 32 -DnaK system and blue for HrcA-GroE system. Fq means the initial value and AF means fold 
change in F(T). Capitals at the bottom (top) of each of the bars show what conditions (A,B,C) 
where broken by further decrease (increase) in parameters. (A) - peak amplitude, (B) - time, and 
(C) - recovering conditions in chaperon production rate (see text for detail). 



Discussion 

We have quantitatively investigated similarities and differences in two heat shock systems: 
one characteristic to gram negative bacteria (e.g. E. coli) and another to gram positive (e.g. 
L. lactis). Remarkably, although the two are very different at the level of promoter regu- 
lation, a striking similarity appears at the level of regulatory networks. Both are governed 
by chaperon-mediated negative feedback loops and in both cases chaperon sequestration is 
employed as a stress sensing mechanism. Furthermore the similarity continues at the level 
of the response dynamics - both systems have characteristic rapid transient responses. 

There are three core features characteristic to both systems, that are necessary to generate 
a rapid transient response observed in both systems upon direct heat shock ( a sudden 
increase in unfolded proteins). 

a) The initial rapid increase is governed by chaperone independent rates, which are the 
a 32 synthesis rate or the rate of HrcA conversion into inactive form. The initial slope 
of increase in chaperons is governed by respectively rj and j c . 
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b) The rapid recovery to steady state is governed by chaperone mediated processes (degra- 
dation of a 32 or activation of HrcA) 

c) The peak is the result of a two rather different time scales involved: a rapid dynamics 
of TF (determining rapid increase and decrease) and a slow chaperone synthesis, de- 
termining the time of the peak in transient response, i.e. time when there is enough 
chaperones to deal with increased amounts of unfolded proteins. 

Furthermore, different realizations of negative feedbacks - one through an transcrip- 
tional activator another through transcriptional inhibitor - impose distinct constrains on 
chaperone-TF binding affinities. Our analyses predicts that whereas the tighter TF- 
chaperone binding increases dynamic range for a 32 system, it would work in opposite di- 
rection and decrease the dynamic range for HrcA system. In other words, chaperone-TF 
binding affinity imposes a lower limit on the amounts of free chaperones for a 32 system, 
where it becomes the upper limit for HrcA system. With the experimentally determined 
binding affinity for a 32 system, the lower limit for free chaperons is 1 nM (see figure H|), 
i.e. one free chaperone per cell, which is low enough to account for possible variations in 
chaperone levels. (The upper bound in this case is determined by the amounts of TF to be 
not less than 1 — 2 protein per cell (1 — 2 nM) so that each cell feels change in TF, thus 
setting upper limit to about 10 4 , see figure HI) 

We predict that the chaperone-HrcA binding should match the upper limit of the desired 
amount of free chaperons and when measured can thus serve as indirect indication of the 
amounts of free chaperones in the cell. 

Acknowledgement 

We thank Kim Sneppen and Steen Pedersen for useful discussions. The research was 
supported by the Danish National Research Foundation. MI is supported by JSPS Research 
Fellowships for Young Scientists and AT us supported by Steno fellowship. 



[1] Dietmar Kltz. Evolution of the cellular stress proteome: from monophyletic origin to ubiqui- 
tous function. J Exp Biol, 206(Pt 18):3119-3124, Sep 2003. 

15 



[2] E. A. Permina and M. S. Gelfand. Heat shock (sigma32 and hrca/circe) regulons in beta-, 
gamma- and epsilon-proteobacteria. J Mol Microbiol Biotechnol, 6(3-4): 174-181, 2003. 

[3] F. Narberhaus. Negative regulation of bacterial heat shock genes. Mol Microbiol, 31(l):l-8, 
Jan 1999. 

[4] H. El-Samad, H. Kurata, J. C. Doyle, C. A. Gross, and M. Khammash. Surviving heat shock: 
control strategies for robustness and performance. Proc Natl Acad Sci USA, 102(8) :2736- 
2741, Feb 2005. 

[5] K. Bourke Arnvig, S. Pedersen, and K. Sneppen. Thermodynamics of heat-shock response. 

Phys Rev Lett, 84(13):3005-3008, Mar 2000. 
[6] M. Wetzstein, U. Vlker, J. Dedio, S. Lbau, U. Zuber, M. Schiesswohl, C. Herget, M. Hecker, 

and W. Schumann. Cloning, sequencing, and molecular analysis of the dnak locus from bacillus 

subtilis. J Bacteriol, 174(10) :3300-3310, May 1992. 
[7] J Gamer, G Multhaup, T Tomoyasu, JS McCarty, S Rudiger, HJ Schonfeld, C Schirra, H Bu- 

jard, and B Bukau. A cycle of binding and release of the dnak, dnaj and grpe chaperones 

regulates activity of the escherichia coli heat shock transcription factor sigma(32). Embo 

Journal, 15(3):607-617, February 1996. 
[8] D. B. Straus, W. A. Walter, and C. A. Gross. The heat shock response of e. coli is regulated 

by changes in the concentration of sigma 32. Nature, 329(6137) :348-351, 1987. 
[9] Steen Pedersen, Phillip L. Bloch, Solvejg Reeh, and Frederick C. Neidhardt. Patterns of 

protein synthesis in e. coli: a catalog of the amount of 140 individual proteins at different 

growth rates, May 1978. 

[10] Niels Ole Kjeldgaard Ole Maale. Control of Macromolecular synthesis. W.A Benjamin Inc, 
New York, 1966. 

[11] W Schumann. The bacillus subtilis heat shock stimulon. Cell Stress & Chaperones, 8(3):207- 
217, 2003. 

[12] C. Lopez-Larrea. Sensing in Nature. Advances in Experimental Medicine and Biology. 
Springer, 2012. 

[13] Jitka Holubov and Jytte Josephsen. Potential of abis as defence mechanism determined by 
conductivity measurement. J Appl Microbiol, 103(6):2382-2391, Dec 2007. 

[14] T. Wiegert and W. Schumann. Ssra-mediated tagging in bacillus subtilis. J Bacteriol, 
183(13) :3885-3889, Jul 2001. 

16 



[15] Eric Guisbert, Takashi Yura, Virgil A Rhodius, and Carol A Gross. Convergence of molecular, 

modeling, and systems approaches for an understanding of the escherichia coli heat shock 

response. Microbiol Mol Biol Rev, 72(3):545-554, Sep 2008. 
[16] In the numerical simulation of the model, the complex formation were solved by using ordinary 

differential equations with much faster time scales than eqs. jl]) to ([3]) for the simplicity of 

calculation. 

[17] From eqs. ([T]) and (jH), the degradation term is given as jci&t]/ l[D /])■ Noting Kj/[Df] ~ 
1/10 (see Table |T] and figure H|) , l/"f c approximately gives the lifetime of <r 32 . 

[18] In the original experiments by Arnvig et al., the response to temperature shift from 30 to 
37°C was peaking at about 0.1 generations, corresponding to 3 — 4 minutes. To address other 
temperature shifts, the time in our simulations can be rescaled in terms of cell generations 
such that the key criteria and the main results would still hold. 



17 



